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Abstract 

This paper is the second of two papers devoted to the study of the evolution of the cosmological horizons 
(particle and event horizons). Specifically, in this paper we consider the extremely general case of an accel- 
erated universe with countably infinitely many constant state equations, and we obtain simple expressions 
in terms of their respective recession velocities that generalize the previous results for one and two state 
equations. We also provide a qualitative study of the values of the horizons and their velocities at the 
origin of the universe and at the far future, and we prove that these values only depend on one dominant 
state equation. Finally, we compare both horizons and determine when one is larger that the other. 

Keywords: Cosmological Parameters - Cosmology; theory - Cosmology: dark energy 



This is an author-created, un-copyedited version of an article accepted for publication in 
Journal of Cosmology and Astroparticle Physics. lOP Publishing Ltd/SISSA Medialab srl 
is not responsible for any errors or omissions in this version of the manuscript or any version 
derived from it. The definitive publisher authenticated version is available online at: 

http : //iopscience . iop . org/1475-7516/2013/02/015 



1 Introduction 

For any observer A we can define two regions in the instantaneous three-dimensional space t = to. The 
first one is the region defined by the events that have already been observed by A, and the second one is its 
complement in the three-dimensional space, i.e. the region that cannot be observed by ^4 at a time to. The 
boundary between these two regions is the particle horizon at to, that defines the observable universe 
for A. Notice that the particle horizon takes into account only the past events with respect to A. Another 
horizon could be defined taking into account also the A's future. This horizon is the event horizon and 
it is defined as the hyper-surface in space-time which divides all events into two classes, those that will be 
observable by A, and those that are forever outside A's range of observation. This horizon determines a 
limit in the future observable universe [6]. 

A deep study of the horizons has been made before in [2] where only one state equation was considered. 
However, it is currently widely accepted the concordance Universe, i.e. that our present Universe is dom- 
inated by two state equations (cosmological constant and dust), one of them of negative pressure, that 
drives the universe into an accelerated expansion. In [4] we obtained a generalization of [2] for the case 
of a concordance universe. This generalization is physically acceptable for all times apart from those near 
to the Big Bang, where the radiation state equation is dominant over cosmological constant and matter 



equations. In this paper we are going to generalize all these results studying the most general case where 
countably many constant state equations are considered in an expanding universe/ Notice that we will 
deal with completely arbitrary state equations as long as the proportion factor remains constant (details 
are provided in the following sections). This condition is in fact not very restrictive, as only the phantom 
energy and some quintessences [1] are accepted to have non-constant state equations. Notice also that the 
particular case with a finite number of state equations is completely covered by our study. 

This paper could be summarized as follows: in section two, we briefly introduce some basic cosmological 
concepts and establish the nomenclature to be used throughout the paper. In sections three and four, 
we define the formal concepts of particle and event horizons, and derive their integral expressions at any 
cosmological time. In section five we obtain the evolution of the horizons in a general universe with 
countably many constant state equations. In section six we compute the values of both horizons at the 
beginning of the universe and at the far future, we also study when one is larger than the other. Last 
section is devoted to gathering all the results obtained, discussion and conclusions. Throughout this paper, 
we use geometrical units where c = G = 1. 



2 Cosmological Prerequisites 
2.1 Proper Distance 

In order to introduce the particle and event horizons, we need to define the proper distance Dp, which is 
the distance between two simultaneous events at a cosmological time to measured by an inertial observer. 
Considering a homogenous and isotropic universe, we can write the Robertson- Walker metric, where a{t) 
is the scale factor (albeit with units of longitude in our system) and k the sign of the curvature: 

ds^ = -dt^ + a{tf ( —^—dr'^ + r-^dQA (2.1) 
\ 1 — fcr^ J 

The definition of the proper distance and some physical considerations [4, for more details see 2.1] allow 
us to express the proper distance as a distance measured using light: 

Dp-'^o ^ (2.2) 
Jt a(i ) 

So to speak, the previous formula represents the distance covered by the light between two points of the 
space-time, but considering an expanding universe since the a factor accounts for this expansion. 



2.2 Hubble parameter 

The Hubble parameter H is defined as: 



i/=i^ (2.3) 

a dt ^ ' 

whose value at t = to is the Hubble constant Ho- Notice in particular that for an expanding universe, a 
is a strictly increasing function, then H never vanishes. The product of H times the proper distance, has 
dimensions of velocity and is known as the recession velocity: 

Vr = HDp (2.4) 

which physically is the instantaneous velocity of an object at a distance Dp with respect to an inertial 
observer. 

We are considering the universe as a perfect fiuid with density p and pressure p, but it can also be 
approximated as composed by non- interacting constituents of density pi and partial pressure pi, where 
all together add up to p and p respectively. Each of this constituents is ruled by its own state equation 
Pi = Wipi. Using the Friedman equations obtained from the Einstein Field Equation and the Robertson- 
Walker metric, and considering the i-th state equation only, we can obtain [3, chap. 3] the i-th density in 
terms of the redshift z: 

p,=Pio{l + zf'+^^^ (2.5) 

Now we define the following magnitudes: 



^Let us recall that countable stands for both finitely many or countably infinitely many, though in our case the finite case 
trivially included in the infinite one. 
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• Critical density as pc = —H^ 

• Dimensionless energy density as Q. = — 

pc 

pi 

• Dimensionless i-th energy density as f2i = — 

pc 

All this quantities, when referred to the current time, are denoted with a zero subindex. If we now substitute 
equation (2.5), pc and Q,iQ in the definition of Q.i , we obtain: 

Finally, adding all the fii (where the index i ranges over all possible state equations) leads to: 



H{z) = Ho 



oo 

xElf(l + -^ (2.7) 

where rii = 3(1 + Wi) are real numbers that can depend on cosmic time or, equivalently, on z. However, 
in this work we will consider that they are all constant, which is indeed a very good approximation for 
the most commonly considered state equations: non-relativistic matter {w = 0), radiation (w = 1/3), 
cosmological constant (lo = —1), or some quintessences (—1 < w < —1/3). 



Remarks. 

2.1 It is important to notice that the currently accepted theories admit only —1 < uii < 1 [1], although 
in what follows we will not use this fact and we will consider the most general situation where Ui are 
arbitrary constant. 

2.2 The case with finitely many state equation is included in the previous expression just by taking 
Q.iO = for every i > K, with K the number of state equations. 

2.3 We could define the density of curvature as Q.c = 1 — and include it as a constant state equation. 
Then, when adding all the Q.i (including Sic), they all add up to 1. Then we can take without loss of 
generality Q. = 1. 

2.4 We will assume the following hypothesis on H (where Sn {z) is the finite summation up to the term 
n and F the infinite summation): 

The derivatives of Sn converge uniformly to a function G on / = [— 1, cxd), that means that: 

sup \S'„{z) - G{z)\ > 

2.5 We can easily weaken this condition by asking the uniform convergence on [—1, b] for every 6 > — 1 
(notice that this does not imply the uniform convergence on [—1, cjo)), but we have to be careful when 
defining the function G. 

2.6 This hypothesis, together with the fact that the series converges pointwise at 2: = 0, implies two 
important facts. First, it assures us that the function H is well defined on [—1, cxj), and in the second 
place it ensures its differentiability and that its derivative can be obtained differentiating term by 
term [5, 9.41]: 

OG 

F'(2) = ^n,n«(l + 2)"'^' 

i = l 

2.7 Notice that this hypothesis is obviously attained if we consider finitely many state equations. 



3 Particle horizon 

3.1 At the present cosmic time 

As the age of the Universe and the light velocity have finite values, there exists a particle horizon Hp, that 
represents the longest distance from which we can retrieve information from the past, so it defines the past 
observable universe. Then the particle horizon Hp at the current moment is given by the proper distance 
Dp covered by the light coming from t = Q (the origin of the universe^) to Iq: 

■^Where a hot big bang is assumed, compatible with the dominant equation of state. 
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Hp = lim Dr, (3.1) 



It is a well known fact that the scale factor a is related with the redshift z as: 



1 + 2=^ (3.2) 

Applying the previous equation, its derivative, the definitions of H and Hp, it can be proved [4, for more 
details see 3.1] that 

'^y (3.3) 



Ho Jo VEf^.o(l + y)"- 
3.2 At any cosmic time 

We have so far obtained the particle horizon for the current cosmic time (corresponding to 2 = 0). If 
we want to know the expressions of this horizon at any cosmic time, or equivalently at any other z' , we 
have to replace in (3.3) the constant parameters Ho and flio (corresponding to 2: = 0) with the parameters 
corresponding to z' that we denote H{z') and Qi{z'). Where we recall from equation (2.6) that ^li{z') is 
i-th dimensionless energy density given by: 

^^{z') = n,oj0^{l + z'r (3.4) 

where Ui — 3(1 + Wi). Therefore, if we replace all these constants in equation (3.3) and simplify the result, 
we obtain the particle horizon for any redshift z (where the prime is omitted in order to simplify the 
notation): 

Hp{z) = ^ r '^^ (3.5) 

Ho Jo VEfi.o(l + ^)"'(l + t/)"^ 

That expression can be rewritten in a more suitable way by making the change of variable given by 
a; = (1 + z){l + y) — 1, then considering the definition of H given by eq. (2.7) we obtain: 

To obtain a characterization of the convergence of this last integral we need the following result of basic 
calculus [5, sect. 7.1]: 

Theorem 3.1. 

Let I = [a, b) where b > a and possibly b — 00, and let /, g be positive and continuous functions in I such 
that 

lim 44=i>0 

x^b- g(x) 

then f{x)dx and g{x)dx either both converge or both diverge. 

Obviously, we have analogous results for / = (a, b]. In what follows, we will denote A'^ the highest Ui which 
appears in H{z) and m the lowest one (we allow them to be ±00 if the set {nijigN is not upper/lower 
bounded). 

If we take f{x) = H{x) and g{x) = (1 + x)'^^^, then its ratio tends to Ho\/Tlm) > when x — >■ 00, so the 
integrals of 1 over these functions behave in the same way. It can be easily proven ]5] that for any 2 > — 1, 
the p- integral of 1/g converges if and only if N/2 > 1. Then we have the following characterization: 



Hp{z) exists for 2; > — 1 if and only if TV > 2 



The limit cases z = —1 and z = ca are studied in section 6. 



4 Event horizon 

4.1 At the present cosmic time 

The event horizon represents the barrier between the future events that can be observed, and those that 
cannot. It sets up a limit in the future observable universe, since in the future the observer will be able to 



obtain information only from events which happen inside their event horizon. According to its definition 
the event horizon can be expressed as: 



He = hm {-Dp) 



(4.1) 



Proceeding as in section 3.1 [4, see 4.1 for more details] we obtain that the event horizon 



dy 



Ho J-i + 

4.2 At any cosmic time 

Analogously as in section 3.2, we can obtain that the event horizon at any z is: 

1 f~ dx 



(4.2) 



H,{z) 



l + z 



(4.3) 



^1 Hix) 

If we now take f{x) = H{x) and g{x) = (1 + a;)™''^, then its ratio tends to HqV^no > when x — >■ — 1, 
so the integrals of 1 over these functions behave in the same way. Now, applying the p-convergence [5], we 
have that for any z > —1, the integral of 1/g converges if and only if m/2 < 1. Then now we have the 
following characterization: 



Hs(z) exists for z > — 1 if and only if m < 2 



The limit cases z = —1 and z = ca are studied in section 6. 
Remark 4.1. 

Notice that if we consider just one state equation, then N = m. Therefore we recover the well known 
result that there cannot exist both horizons at the same time, and in fact, if A'^ = ?n = 2 (state equation 
of curvature or cosmic strings) there exists none of them. 



5 Study of the evolution of the horizons in a general universe 

In this section we obtain a general formula for the evolution of the horizons whenever they are defined. In 
order to do that we just need the following results: 



F{x) 
dz 

IE ^ 



b(x) 



f{t)dt ^ ^ = f{b{x))h\x) 



Ho da ao n , \tj-/ \ 



(5.1) 
(5.2) 



Where the first expression follows from the well known fundamental theorem of calculus [5, sect. 7.1] where 
/ is a continuous function defined in a close interval, and the last one is obtained applying first equation 
(3.2), then equation (2.3) and finally equation (3.2) again. Now deriving equations (3.6) and (4.3): 



dHp 
dz 



dHe 

dz 



(1 + ^)2 

Hp 

l + z ^ 
1 

ll + zf 



dx 1 
H{x) ~ {l + z)H(z) 
1 



[l + z)H{z) 
dx 



+ 



1 H{x) {l + z)H{z) 
1 



1 + 2 ■ {l + z)H{z) 
Now, taking into account the chain rule and equation (5.2), we obtain: 



dH , 
dt 



^ = Hp{z)H{z) + 1 



dHe 
dt 



= He{z)H{z) - 1 



(5.3) 



Note that Hp{z)H{z) and He{z)H{z) represent, respectively, the recession velocities of the particle and 
event horizons (2.4). Its physical meaning is explained in subsection 7.2, where all the conclusions are 
provided. 
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6 Relevant values of the horizons and their velocities 



6.1 Horizons at the Origin 

In this section we make a deep study of the behaviour of the particle and event horizons, as well as its 
derivatives at the origin of the universe z = <x. From now on, when dealing with the particle horizon, we 
may assume that it is defined for 2 > — 1 i.e. that N > 2. Similarly, when dealing with the event horizon 
we may assume that m < 2. 

As Hp exists, then the integral on equation (3.6) is convergent for any z > —1, so if we take the limit 
— >■ 00, Hp goes to zero. Then 

lim Hp(z) = (6.1) 

2—^00 

If we try to compute the limit of He{z) when z — cxj, we obtain zero if A'^ > 2 (as in this case, the 
integral appearing in the expression of He is convergent), and an indetermination ^ if A < 2. For 
the latter case, applying L'Hopital's rule and the fundamental theorem of calculus, we obtain: 



Hp 



He 



lim HAz) = lim — 

c— )-oo z—^oa M i^Z 



if AT > 

1 



00 if A < 



if A = (6.2) 



Once we have computed the limits of the horizons themselves, let us deal with their velocities i.e. how the 
proper distance of the horizon varies with respect to the time coordinate t. In order to do that, first notice 
that: 



if AT > 2 



lim ^ = I HoV^i if TV = 2 (6.3) 
if TV < 2 

^'W = ^E^>"«(1 + ^)"'"' (6.4) 
Where the last formula is valid thanks to the hypothesis we assume on H (see remarks after eq. (2.7)). 
As we are assuming A" > 2, if we take the limit 2; — > 00 of HpH we obtain an indetermination • co: 

dx L'H 



lim HpH = lim — — tjtt^ . , „, . 

z->oo z^oo (I ^ z)/ H{Z) J ^ H(x) 



lim „ ~ lim ^ 



z^oo H-(l + z)H' z^oa H ~ (1 + Z)H' 

lim ^1±^ - 1^ 



1 ,. J] ?iif2io(l + z) 
- lim 



2z^^ Y,Q.m{\ + zY^ 



Then, from equation (5.3) we have 



2 7 A- 



(6-5) 



We are now assuming that m < 2, but when z tends to infinity, the dominant term is the one with 
the largest Ui i.e. the one associated to A, and so we have to consider three possibilities: 
A > 2 In that case, 

H{x) 

and equation (6.3) implies that velocity of the event horizon tends to infinite. 



- 6 - 



N = 2 Now the previous limit is divergent, but again equation (6.3) implies that velocity of the event 
horizon tends to infinite. 

N < 2 Again, we have that the previous integral is divergent, but in this case equations (5.3) and (6.3) 
lead to an indetermination • cxd. Using exactly the same computations as in the dHp/dt case 
(but now, when applying the fundamental theorem of Calculus no minus sign appears) we have: 

1-^ = Af (6-6) 

dt 2- N ^ ' 

In table 1 we have summarized all these data. 

6.2 Horizons at the far future 

In this section we make a deep study of the behaviour of the particle and event horizons, as well as its 
derivatives at the far future z = —1. 



Hr, 



Notice that the integrand on equation (3.6) is positive everywhere (see explanation below eq. (2.3)), 
and can only vanish for the value x = z — —1, then the integral exists (though possible infinite) and 
it is greater than zero, so the fraction multiplying the integral in (3.6) leads to 

lim Hpiz) = oo (6.7) 

2— > — 1 

If we try to compute the limit of He, and indetermination 2 appears. Proceeding in analogy with the 
previous section, we obtain: 

fO if ?n< 
^= if m = (6.8) 
oo if m > 
Now, in order to compute the limits of both velocities, notice that 

oo if m < 2 

lim = J HoV^> if m = 2 (6.9) 



z-*-! I + Z 



if m > 2 



dHp 



dt 



We are assuming N > 2, but when z tends to —1, the dominant term is the lowest one. Making the 
same study as in the dHe/dt case for z oo we obtain: 

(oo if m < 2 
m 
if m > 2 
m — 2 

We are now assuming that m < 2, and if we now proceed as in the dHp/dt case for z — >■ oo (as an 
00 • indetermination arises) we have: 

lim ^ = (6.10) 

t->cx) dt 2 — m 

In table 1 we have summarized all these data. 
Remarks. 

6.1 Whenever it exists. Hp always comes from zero at the origin of the universe, and tends to infinite at 
the far future no matter whatsoever which and how many state equations are considered. 

6.2 Notice that He is constant for every z G ( — 1, oo) if and only if just one state equation with N = m = 
is considered, in which case, He{z) — Ho^/Tlo for every z. 

6.3 Though physically it is not accepted in the current theories, we can deduce that if m < or equiv- 
alently Wm < —1, we would have that the event horizon speed becomes strictly negative for some 
z > — 1 as in the limit it is strictly negative. We proved in [4] that for the concordance Universe (two 
state equations with m = and A'' = 3) such a thing does not happen i.e. it always remains positive. 

6.4 Notice that all the results we have obtained so far, depend on the state equations associated to A*" or 
m, and some of them do not depend on the state equations considered. 
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(m < 2) 



dH, 



dt 



^ (TV > 2) 



dH, 
dt 



(to < 2) 



Origin of 
the universe 

t = 



Far Future 

t ~ oo 







1 



Hoy/ ^ NO 

oo 







1 



HoV^mO 
CO 



if > 
if A^ = 
if < 
if TO < 
if ?71 = 
if TO > 



N 



N -2 



oo if m < 2 

TO 



TO- 2 



if TO > 2 



oo if iV > 2 

N 



2- N 



if iV < 2 



2 — TO 



Table 1 - Some important values of the horizons and their velocities in geometrical units. 



6.3 When are Hp and H^. equal? 

It is interesting to know whether both horizons might have or not the same values, and if so, how often this 
could happen. In order to do that, let us look for a zo such that Hp{zo) = -ffe(zo), where we are obviously 
assuming that N > 2 and m < 2 so both horizons exist. First notice that if zq = —1, then according to 
equation (6.7), they have both to be infinite and in that case we know from equation (6.8) that m > 0. 

The previous solution has little interest, so let us now look for other possible solutions such that zq > —1. 
According to equations (3.6) and (4.3), and cancelling out the fractions before the integrals: 



r°° dx _ r 



"0 dx 



Hix) J_, H{x) 
The fact that N > 2 and m < 2 ensures that the integral from z 

dx dx 



L ; 



1 H{x) 



1 H{x) 



+ 



— 1 to infinite is convergent, then: 
dx 



H{x) 



= {l + zo)[H^{zQ) + Hp{zQ)'j = 
(l + zo)2He(^o) " 



/-I H{x) 

Where in the ★ equality we have used that Hp{zo) — He{zo). Then finally: 



i: 



dx 



L 



As H{x) is positive wherever it is defined, the integral of the left hand side of the equation is strictly 
increasing with respect to zq, and it goes from when zo — —1, and tends to L when zo ~> cxd. Then 
there exists a, zo G (— l,oo) such that verifies the previous equation (6.11) and then verifies the equation 
Hp{zo) = Hs{zo)- If we suppose that there exist zi and Z2 verifying the previous equation, then: 



dx _ r 



dx 
dx 



1 H{x) H{x 



'-■^ dx 
1 H{x) 

1 dx 



2 2 



and again, as H is positive, we have that zi = Z2. Finally, taking into account the following chain of 
equivalences: 



Heiz) > Hp{z) 



i 



dx 

.mx) 

^ dx 



dx 



dx 

-1 m^) 



dx 

" dx 



'° dx 
> ^ 



Z > Zo 
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where in the last equivalence we have use once again that H is positive, we may summarize all the previous 
computations in the following result: 

Theorem 6.5. 

There exists one and only one zo > —1 such that 

He{z) > Hp{z) for every z £ {zo,oo) 
He{zo) = Hp{zo) 

He{z) < Hp{z) for every z e {—l,zo) 

Remark 6.6. 

Obviously, the numerical value of the zo itself depends on which state equations we take into account and 
also on the numerical value of all the Qio- 



7 Results and interpretation 

7.1 Application to the concordance universe 

In this section we are going to prove that in fact we are generalizing the previous results obtained for a 
concordance universe in [4] i.e. taking into account two state equations (dust and cosmological constant). 
For these two state equations we have that A'' = 3 and m = 0, then table 1, recovering the constant c, 
becomes table 2: 





z t 






dHp 
dt 


dHe 
dt 


Origin of the universe 


oo 








3c 


oo 


Future 


— 1 oo 


oo 


c 


oo 






Table 2 - Relevant values for the concordance universe. 



which agrees completely with table 1 of [4], where we considered the concordance universe i.e. euclidean, 
with cosmological constant and dust. 

We can go further in our analysis as it is widely accepted that for the origin of the universe, the state 
equation that dominated was radiation (then we have = 4 instead of A'^ = 3). In this case it holds that 
all the data are the same as for the concordance universe, but for the velocity of the particle horizon at 
the origin (see table 3): 

Um = 2c 

t^o dt 





z t 


Hp{z) 


Heiz) 


dHp 
dt 


dHe 
dt 


Origin of the universe 


oo 








2c 


oo 


Future 


— 1 oo 


oo 


c 


oo 





Hoy/rimO 



Table 3 - Relevant values for the concordance universe with radiation. 



7.2 Conclusions 

In this paper, we have first obtained integral expressions for the particle and the event horizons considering 
the extremely general case of a universe with or without curvature ruled by countably infinitely many state 
equations (under some reasonable hypotheses on those state equations that are fulfilled when only finitely 
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many state equations, no matter which one, are considered). After some suitable manipulation, we have 
computed their first derivative, obtaining these extremely simple equations; 

^ = H,{z)H{z) + 1 ^= H4z)H{z) ~ 1 

where HpH and HeH are the recession velocities of the particle and event horizons respectively. These 
expressions generalize the previous stated results derived for just one constant state equation and null 
curvature [2] and for two state equation [4]. Notice that all these results and values are independent of the 
current values of the involved constants. 

The equations of the velocity of the horizon have a remarkable physical meaning. As the recession velocity 
HpH is the instantaneous velocity of an object located at the distance of the particle horizon Hp, from 
the first equation we deduce that the instantaneous velocity of the surface of the horizon particle is faster 
than the one of the objects over the particle horizon, then more and more objects are entering into the 
particle horizon and they will never get out. Analogously, the second equation stands that the recession 
speed of the event horizon is slower than that of the objects over the event horizon, and then more and 
more objects are disappearing from the event horizon and they will never get into again. 

We have also obtained the value of the horizons and their velocities at origin of the universe and for the 
far future with the great generality stated above. As showed in table 1, some of them depend on which 
state equations are considered. We have obtained also the expected result, that these values only depend 
on the respective dominant state equation. Finally we have obtained that there exists one and only one 
20 > —1 where both horizons have the same value, and also that for lower values of the redshift z it holds 
He{z) < Hp{z) while for larger values He{z) > Hp{z). 
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